Anisotropic Magnification Distortion of the 3D 
Galaxy Correlation: I. Real Space 
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It has long been known that gravitational lensing, primarily via magnification bias, modifies the 
observed galaxy (or quasar) clustering. Such discussions have largely focused on the 2D angular 
correlation function. Here and in a companion paper (Paper II) we explore how magnification bias 
distorts the 3D correlation function and power spectrum, as first considered by Matsubara (2000). 
The interesting point is: the distortion is anisotropic. Magnification bias in general preferentially 
enhances the observed correlation in the line-of-sight (LOS) orientation, especially on large scales. 
For instance, at a LOS separation of ~ 100 Mpc/h, where the intrinsic galaxy-galaxy correlation 
is rather weak, the observed correlation can be enhanced by lensing by a factor of a few, even at 
a modest redshift of z ~ 0.35. This effect presents an interesting opportunity as well as challenge. 
The opportunity: this lensing anisotropy is distinctive, making it possible to separately measure the 
galaxy- galaxy, galaxy-magnification and magnification-magnification correlations, without measur- 
ing galaxy shapes. The anisotropy is distinguishable from the well known distortion due to peculiar 
motions, as will be discussed in Paper II. The challenge: the magnification distortion of the galaxy 
correlation must be accounted for in interpreting data as precision improves. For instance, the 
~ 100 Mpc/h baryon acoustic oscillation scale in the correlation function is shifted by up to ~ 3% 
in the LOS orientation, and up to ~ 0.6% in the monopole, depending on the galaxy bias, redshift 
and number count slope. The corresponding shifts in the inferred Hubble parameter and angular 
diameter distance, if ignored, could significantly bias measurements of the dark energy equation 
of state. Lastly, magnification distortion offers a plausible explanation for the well known excess 
correlations seen in pencil beam surveys. 

PACS numbers: 98.80.-k; 98.80.Es; 98.65.Dx; 95.35.+d 



I. INTRODUCTION 

Since the pioneering work of 0, it has been appreci- 
ated that gravitational lensing modifies the spatial dis- 
tribution of high redshift objects, such as galaxies or 
quasars (henceforth, the term 'galaxies' can be viewed 
as synonymous with quasars or any sample of survey ob- 
jects). Suppose for instance there is a large mass con- 
centration between some galaxies and the observer. The 
observed number density of galaxies would decrease due 
to the stretching of the apparent inter-galaxy spacing, 
but would increase due to the enhanced ability to see 
very faint galaxies that otherwise would have gone unde- 
tected. The size of the net effect depends on the number 
count slope. This effect is known as magnification bias. 

The implications of this effect for the observed galaxy 
angular correlation function were worked out in 0, [H, 
0, HI- It was pointed out by [(| that the magnifica- 
tion bias correction can be isolated by measuring the 
angular cross-correlation between galaxies at widely sep- 
arated redshifts. This idea was subsequently realized in 
measurements by @, H| through the cross-correlation of 
quasars and galaxies. See [|| Gil [III f° r some recent theo- 
retical work and references therein for earlier discussions 
of quasar-galaxy associations. See also [12] on the early 
use of magnification bias to make mass maps of galaxy 



clusters. 

With the important exception of Matsubara (2000) [l3| 
(to which we will return at the end of this section), earlier 
papers have focused on the 2D angular correlation func- 
tion. Here and in Paper II [l4[ of this series, we study the 
effect of magnification bias on the 3D correlation function 
and power spectrum, ft is not hard to convince oneself 
that lensing makes the 3D correlation anisotropic. There 
are several effects at work, but the simplest one to think 
about is the following. Correlation function is measured 
by pair counts of galaxies. A pair of galaxies that are 
aligned along the line-of-sight (LOS) behave differently 
from a pair oriented transverse to the LOS. In the for- 
mer case, the closer galaxy can lens the background one. 
The same does not happen in the transverse orientation. 
The net effect is an anisotropic observed correlation. One 
might think such an effect must be small: after all, typi- 
cal LOS separations in clustering measurements are much 
smaller than the depth of surveys i.e. the lens is located 
much closer to the source than to the observer. However, 
one must remember that the lensing effect grows with the 
LOS separation, whereas the intrinsic galaxy correlation 
generally drops with separation. At a sufficiently large 
LOS separation, e.g. ~ 100 Mpc/h, where the intrinsic 
galaxy correlation is rather weak, one should not be too 
surprised that the lensing induced correlation can actu- 



2 



ally dominate. As we will show, the magnification bias 
induced anisotropy has a distinctive shape: 

Us{Sx,8x±)= i gg {^J 5 X 2 + 8x\) (1) 
+f(Sx±)Sx + g(Sx±) 

where 5\ and Sx± are the LOS and transverse separa- 
tions respectively, £ b s is the observed correlation, is 
the intrinsic galaxy-galaxy correlation, fSx is the galaxy- 
magnification correlation and g is the magnification- 
magnification correlation. Note that / and g are func- 
tions of the transverse separation only. 

This distinctive anisotropy pattern makes it in prin- 
ciple possible to completely separate the three differ- 
ent contributions from data: galaxy-galaxy, galaxy- 
magnification and magnification-magnification correla- 
tions. (We will discuss other sources of anisotropy, such 
as peculiar motions in Paper II.) In a sense, the earlier 
work @, 0, H, OH [H| on angular correlation func- 
tion for galaxies (or galaxies and quasars) at widely 
separated redshifts focused on one particular limit: a 
very large LOS separation Sx such that the galaxy- 
magnification correlation fSx dominates. Studying more 
moderate LOS separations allows one to measure the 
magnification-magnification correlation g, which is per- 
haps of more theoretical interest since it relates directly 
to the mass. The key is to use the full 3D information, 
i.e. exploit the distinct dependence on Sx and Sx± of 
each term in eq. (JTJ to separately determine all three 
correlations. 

Another important implication of eq. ([T]) is that 
care must be taken in interpreting galaxy clustering 
data. For instance, future galaxy surveys hope to de- 
termine the baryon oscillation scale to high precision 
ES US El HI, 01 M, HI ■ Lensing induced shifts of 
the apparent baryon oscillation scale in the LOS direc- 
tion affects the inference on the Hubble parameter H(z), 
while shifts in the transverse direction affects the infer- 
ence on the angular diameter distance. Recall that at a 
redshift of 1, a ~ 1% shift in angular diameter distance, 
or a ~ 1% shift in the Hubble parameter, corresponds to 
a ~ 5% shift in the dark energy equation of state. This 
means that even small lensing corrections are in principle 
a worry, not to mention potentially large corrections at 
high redshifts or at particularly susceptible orientations, 
such as the LOS direction. 

The rest of the paper is organized as follows. In m we 
derive and numerically compute the magnification bias 
distortion, magnification distortion in short, of the ob- 
served galaxy clustering - Hll Al focuses on the overall 
anisotropy of the correlation function while i jll Bl focuses 
on the baryon acoustic oscillations. In £11111 we present 
useful order of magnitude estimates which offer deeper 
insights into the numerical results of [jH] We conclude in 
^IVl with a discussion of longstanding puzzles posed by 
pencil beam surveys. Appendix [A] contains a discussion 
of lensing corrections to the observed correlation other 
than those due to magnification bias, and a generaliza- 



tion of the magnification bias effect to account for more 
complicated galaxy selection. Appendix[B] contains a dis- 
cussion of the higher order Taylor expansion terms that 
are ignored in Hll Al 

In this paper, we focus exclusively on the effects of 
magnification distortion in configuration/real space. Pa- 
per II of this series explores this effect in Fourier space. 
Redshift space distortion due to peculiar motion and the 
Alcock-Paczynski effect will be discussed there as well. In 
yet another paper, we study the effects of magnification 
bias on the angular galaxy correlation function, focusing 
on the impacts on features of the power spectrum [22| . 

As this paper was being completed, two preprints ap- 
peared which discussed some related issues. Vallinotto et 
al. [23[ | explored the impact of lensing, especially magni- 
fication bias, on baryon oscillation measurements. They 
focused on the 2 point correlation function where the 2 
points are at exactly the same redshift. They did not ex- 
amine the full 3D correlation, in particular its anisotropy. 
Their findings are consistent with ours for pair separa- 
tions oriented transverse to the LOS, and are more con- 
nected to our other paper on the angular correlation func- 
tion [22j]. Wagner et al. [H| examined the anisotropy of 
the 3D correlation that is introduced by light cone effects. 

After this paper was initially circulated as a preprint, 
a pioneering paper by Matsubara (2000) [l3[ was kindly 
brought to our attention, where he derived an expression 
for 3D correlation function in the presence of magnifica- 
tion bias as well as redshift and cosmological distortions. 
In this paper, we have extended his analysis in a number 
of ways. By showing the 3D correlation function in terms 
of the LOS and transverse comoving separations (as op- 
posed to redshift and angular separations as in [LI), the 
anisotropy pattern is brought out explicitly. We also em- 
phasize the possibility to completely separate the three 
different contributions to the observed correlation func- 
tion (see Fig. [5] below). As an application, we show how 
magnification distortion impacts baryon acoustic oscilla- 
tion measurements. In Paper II, we also study the ap- 
pearance of magnification distortion in Fourier space (as 
well as redshift space), which turns out to have some im- 
portant and interesting qualitative differences from the 
appearance in configuration/real space. 

II. MAGNIFICATION DISTORTION 

Magnification bias introduces a well-known correction 
to the observed galaxy overdensity [l], H, H, 0, @| : 

Sobs = Sg + (2) 

where the observed galaxy overdensity <5 bs, the origi- 
nal/intrinsic galaxy overdensity 5 g and the magnification 
bias correction <5 M are all functions of the galaxy posi- 
tion, specified by the radial comoving distance x an d the 
angular position 6. The magnification bias correction is 
given by 

S„ = (5s - 2)k (3) 
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where k is the lensing convergence: 



x' 



(4) 



where is the gravitational potential, and is the 2D 
Laplacian in the transverse directions. We assume a flat 
universe - generalization to an open or a closed universe 
is straightforward. The symbol s stands for 



dlog w N(< m) 
dm 



(5) 



where N(< m) is the cumulative number counts for 
galaxies brighter than magnitude m. This assumes the 
galaxy sample is defined by a sharp faint-end cut-off. A 
broader definition of s for a more general galaxy selection 
is given in Appendix |Al 

Defining the galaxy bias b by 5 g = b5, where S is the 
mass overdensity, eq. ^ can be rewritten as 



^obs 

~b~ 



= 5 



5s - 2 



(6) 



The relative importance of the intrinsic clustering and 
the magnification bias correction is therefore controlled 
by, among other things, the sample dependent ratio 
(5s — 2)/b. The observed redshift-dependent luminosity 
function [25| can be used together with the halo model 
to estimate this ratio (see [26| for details). Fig. [T] shows 
this as a function of redshift for five different samples, 
each defined by a different B-band apparent magnitude 
cut-off. This figure should be viewed as an illustration of 
the range of possibilities only. The precise values of s and 
b depend sensitively on details of how the galaxy/quasar 
sample is selected, for instance subject to color cuts and 
so on. Unless otherwise stated, we adopt throughout 
this paper the value (5s — 2)/b = 1 to illustrate the ef- 
fect of magnification bias on clustering measurements. 
For the correlation function or power spectrum, one can 
roughly scale the magnification bias correction we obtain 
by (5s-2)/6 for galaxy redshift < 1.5 and by {5s-2) 2 /b 2 
for redshift > 1.5 (the former is dominated by the galaxy- 
magnification cross-term while the latter is dominated by 
the magnification-magnification term; see S jll A|) . Note 
also that we assume a scale independent (linear) galaxy 
bias b. Nonlinear galaxy bias, as we will see, is important 
in certain situations including, surprisingly, some where 
large scale clustering measurements are involved ( mi A\i . 

The reader might wonder whether there are other lens- 
ing corrections to the observed galaxy overdensity. In- 
deed there are, and they are discussed further in Ap- 
pendix [AJ It suffices to say magnification bias is the 
dominant effect for regimes of practical interest. 

In all illustrative examples below, we employ the fol- 
lowing cosmological parameters: the Hubble constant 
h = 0.7, matter density fi m = 0.27, cosmological con- 
stant S1a = 0.73, baryon density S~2& = 0.046, power spec- 
trum slope n — 0.95 and normalization erg = 0.8. We 
employ the transfer function of [2?| , and the prescription 
of [28] for the nonlinear power spectrum. 
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FIG. 1: The ratio (5s — 2)/6, where s is the number count 
slope and b is the galaxy bias, as a function of redshift z for 
five different galaxy samples. Each sample is defined by a B- 
band (4344 angstroms) apparent magnitude cut-off as shown 
for each curve. The curves should be viewed as illustrative 
rather than definitive: only linear galaxy bias is accounted 
for here, and the precise values for s and b depend on how the 
galaxies/quasars are selected e.g. subject to color cuts and so 



A. The Anisotropic Correlation Function 



The observed two-point correlation function is: 



£obs(Xl,0i;X2,02) = (<5obs(Xl,0l)<5 o bs(X2,02)) 

= £ S9 (xi, 0i;xa, 02) + £ 9/1 (xi, #i;x2, o 2 ) 

+fg/*(X2) #2! Xl, 0l) + W(Xl> 015X2, 6 2 ) 



(7) 



where 

£ 99 (l;2) = («5 9 (l)5 g (2)) 
^(2;1) = (^(2)^(1)) 



f flP (l;2) = (5 fl (l)^(2))(8) 

U(i;2) = (W^(2)> 



with the arguments 1 and 2 as shorthands for xi, #i and 

X2,#2- 

Using the Limber approximation, the galaxy- 
magnification cross-term(s) can be written as: 



U(Xi,0i;X2,0 2 ) = -tf 2 f! m (5s-2) 

(i + ^i) (x2 ~ Xl)xi e(xi <x 2 ) 

X2 

J (2ir) 2 9m{ ' 



(9) 
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where ©(xi < X2) is a step function which equals 1 if 
Xi < X2 and vanishes otherwise, H is the Hubble con- 
stant today, fi m is the matter density today (normalized 
by the critical density), Z\ is the redshift corresponding 
to the comoving distance xij and P gm is the (3D) galaxy- 
mass power spectrum, and kj_ is the transverse Fourier 
wave vector. We have used the Poisson equation to relate 
the gravitational potential <fi to the mass overdensity 5: 



V 2 (t) = 3H^n m (l + z)6/2 



(10) 



Note that the speed of light is set to 1 throughout. 

The magnification-magnification correlation, or mag- 
nification auto-correlation, is 

fwOa, *i; X2, 02) = [^Hln m {bs - 2)] 2 (11) 

mi„.( Xl , X2 ) (xi _ % / )x / (X2-M 

d-X 

Xl X2 



(1 



J\2 



(27T) 



ikx-x'(^l— 03) 



where P mm is the (3D) mass-mass power spectrum: it 
is evaluated at the redshift z' which corresponds to the 
integration variable x' ■ 

The focus in the literature has been on £ Q b s as an 
angular correlation: xx an d X2 are typically integrated 
over some radial selection functions, and they can sig- 
nify either two different redshift bins (i.e. angular cross- 
correlation between galaxies/quasars at two different red- 
shifts 6]), or the same redshift bin (i.e. angular auto- 
correlation function [4|). 

Here, let us take a slightly different perspective: think 
of £ bs as a 3D correlation function [l3| . For a galaxy sur- 
vey with redshift information (either spectroscopic red- 
shifts or high quality photometric redshifts), this would 
be a very natural thing to do. Further, suppose one has a 
galaxy survey, or a subsample thereof, that spans some fi- 
nite redshift range such that the radial separation xi ~ X2 
is always small compared to xi or X2- This is a sensible 
assumption since at sufficiently large separations, galaxy 
evolution becomes important and complicates one's anal- 
ysis. Let x be the mean radial comoving distance to 
these galaxies, and z be the associated mean redshift. 
The galaxy-magnification cross-correlation and the mag- 
nification auto-correlation can be simplified as follows: 



ffl/iOd, 015X2, 82) +£ SM (X2, ^25X1, Ox) = 
^ff 2 fi m (5 S -2)(l + z)|x2-Xi| 



(12) 



d 2 k A 



■p gm (z, kx_y- 



k±-x(0i-0 2 ) 



^(Xx,0x;X2, 2 ) = [-Hln m {hs - 2)f 



(13) 



dx 



(x - x')x' 



(1+^') 2 



r (PkA_ , , , V ikx. x '(0i-0 2 ) 

J (27r) 2 ' ' 
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FIG. 2: An illustration of how and in principle 

can be obtained from the observed £ bs- At a given transverse 
separation 5x± (2 examples are shown), £ b s (solid line) at a 
large LOS separation 8\ IS dominated by the magnification 
corrections 2£ gM +£ M)J which have the form f(Sx±)6x+g(Sx±) 
i.e. it is linear in 8\ (dotted line). The extrapolation of this 
dotted line to S\ = gives g or equivalently £ MM . Its slope 
gives / which can be multiplied by 8\ to obtain 2£ 9M . Finally, 
subtracting the dotted line from £ b s yields £ 99 . 



where we have Taylor expanded xi and X2 around x and 
retained the lowest order contributions. It is useful to 
compare these two expressions to the intrinsic (unlensed) 
galaxy auto-correlation, or galaxy-galaxy correlation: 

£ gg {Xx,8 1 ; X 2,8 2 )= (14) 

WV(xi-x 2 ) 2 + x 2 (0i-0 2 ) 2 ) = 

f p (z Z-Vk txi-xs) 

J (2tt) 3 " [ ' ' 

ignoring for now the issue of redshift distortion, which 
will be addressed in Paper II. Note that xi and X2 refer 
to the points corresponding to Xi> &x and X2, #2- 

The observed correlation function is a sum of all three 
correlations above (eq. [7], [12], [13] and [H]). (A dis- 
cussion of higher order corrections to the latter three 
can be found in Appendix [B]) Viewed in this way, the 
anisotropy of the lensing induced corrections is quite 
striking: £ w (l, 2) +£ gM (2, 1) scales linearly with the line- 
of-sight (LOS) separation X2— Xi (i- e - it increases rather 
than decreases with the separation!), and £ PM is indepen- 
dent of the LOS separation. The intrinsic galaxy auto- 
correlation £g g is isotropic and generally decreases with 
separation. 
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FIG. 3: z = 0.35: (a) Contours of constant £ bs (solid) and £ 99 (dotted), left to right: 0.01 (magenta), 0.002 (cyan), 0.001 
(black; triple contours), 0.0005 (blue), (red), -0.0001 (green). The LOS separation is 5\', the transverse separation is Sx±. 
£ 99 is isotropic but £ b s is not. (b) £ 99 (black dashed), 2£ 9M (sloped red dotted), £,j M (flat red dotted) and £ b s (red solid) for a 
separation vector oriented along the LOS. The inset shows the ratio £ bs/£ 99 for several other orientations (solid); dot-dashed 
line shows the ratio of the respective monopoles. Note Sx 2 — 8\ 2 + 5x\. (c) A zoomed in view of £ 99 , £ bs, 5x 2 £, gg , 5x 2 ^ ^, B 
around the baryon wiggle for a separation vector oriented along the LOS. Note how dangerous it is to use 5x 2 £ b s to locate the 
baryon peak, (d) Lower panel shows the monopole of £ b s ; upper panel shows the difference monopole £ b s — £ 99 for several 
different values of (5s — 2)/6. Unless otherwise stated (as in panel d), (5s — 2)/b — 1 throughout. All £'s are normalized by b 2 . 



We can summarize the distinctive lensing induced 
anisotropy to the observed correlation as follows: 

Us(Sx,Sx ± )= Z ga {yl 8 X 2 + 8x\) (15) 
+f(6x ± )Sx + g(Sx ± ) 

where 5\ and 5x± are the LOS and transverse 
separations respectively, fSx represents the galaxy- 



magnification correlation and g represents the 
magnification-magnification correlation. Here, / 
and g are functions of the transverse separation only, 
and are determined by the galaxy-mass and mass-mass 
power spectra. This distinctive form of the anisotropy 
allows us in principle to separately measure £ gg , f and g, 
from which we can infer the galaxy- galaxy, galaxy-mass 
and mass-mass power spectra. For instance, at any given 
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FIG. 4: Analog of Fig. [3] for z = 1. The contours in (a) are: 0.01 (magenta) 
contours), (red) and -0.001 (green; double contours). 



0.002 (cyan), 0.0008 (black), 0.0005 (blue; triple 



5x±, plotting £ bs as a function of the LOS separation 
Sx would reveal a linear contribution at sufficiently large 
Sx's where £ gg is very small. Its slope tells us / and its 
extrapolation to 5\ = tells us g. Subtracting Sxf + g 
from £ bs then yields £ gg . This is illustrated in Fig. [5] 
Fig. [2] should be viewed as a proof of concept only. The 
optimal method for extracting £ gg , £ 9(L1 and from 
realistic data requires more investigation. 

Let us study this distortion of the correlation function 
by magnification bias in more quantitative detail. Fig. 
[3] - [7] are a series of figures showing several interesting 
aspects of the observed and intrinsic correlation func- 
tions for galaxies at mean redshifts of z = 0.35, 1, 1.5, 2, 3. 
These redshifts are chosen to match roughly current and 



future surveys 21|. Here, as is throughout the paper, 
the correlation functions shown in all figures are implic- 
itly divided by b 2 . 

Panels (a) show contours of constant £ gg (dotted) 
and £ bs (solid) as a function of the separation vector 
<5x = xi — X2. Symmetry dictates that the two point 
correlations are functions of only the LOS projection 
<5x = Ixi — X2I (y-axis) and the transverse separation 
5x± = x\6i — 02 1 (x-axis). Note how certain values of 
£ gg or £ bs map to multiple contours: this is because the 
correlation function is not monotonic, both around the 
baryon wiggle and beyond the zero-crossing scale (see 
panels b). Redshift distortion due to peculiar motion is 
ignored, and therefore £ gg is isotropic. The observed cor- 
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FIG. 5: Analog of Fig. [3]for z = 1.5. The contours in (a) are: 0.01 (magenta), 0.002 (cyan), 0.0008 (black), 0.0005 (blue), 
(red) and -0.00005 (green; double contours). 



relation £ Q b s is anisotropic due to the presence of mag- 
nification bias corrections. These corrections are most 
apparent for pair separations oriented along the LOS. 
In this orientation, the galaxy-magnification correlation 
grows linearly with separation (see eq. |12| ) while the 
intrinsic galaxy-galaxy correlation £ g9 drops with separa- 
tion. At LOS separations 70 Mpc/h or larger, £ gM domi- 
nates over ^ gg . As expected, the magnification distortion 
of the correlation function becomes increasingly signif- 
icant at larger redshifts - the distortion is confined to 
small Sx± , s at low redshifts, but diffuses to larger Sx± , s 
at high redshifts. The overall distortion pattern vaguely 
resembles the well known finger-of-god (FOG) effect due 
to virialized motion, in that the contours of constant £ b s 



are elongated in the LOS direction. The precise shapes of 
the anisotropy, however, are rather different - FOG cor- 
rections do not have this linear dependence on the LOS 
separation that £ 3/J has. Moreover, FOG due to pecu- 
liar motion generally does not extend out to such large 
scales. We will examine the net clustering anisotropy ac- 
counting for both magnification bias and peculiar motion 
in Paper II. It is also worth emphasizing that £ gA , could 
have the opposite sign if s < 0.4, in which case the LOS 
correlation is enhanced in the negative direction. As can 
be seen from Fig. [U unless one has an exceptionally faint 
magnitude cut-off, s is generally larger than 0.4 at high 
redshifts where gravitational lensing is most effective. 
Panels (b) show £ b s and its three different contri- 
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FIG. 6: Analog of Fig. Ofor z = 2. The contours in (a) are: 0.01 (magenta), 0.002 (cyan), 0.0005 (black), 0.0003 (blue; triple 
contours), (red) and -0.00005 (green; double contours). 



butions for a separation vector that is oriented along 
the LOS (0 X = 0, where cos8 x — Sx/Sx with Sx — 
\Sx\ = ^/6 x 2 + Sx' 2 L ). Eq. (H3D tells us £ w (flat dot- 
ted line) is independent of the LOS separation, and is 
therefore simply a constant for this particular orienta- 
tion. As discussed earlier, £ gM (sloped dotted line) grows 
linearly with separation in this orientation. The result- 
ing Cobs = i ag + 2£ gAl + £ w is quite a bit larger than 
£ gg , even at a redshift as low as z — 0.35. At a scale 
of 100 Mpc/h, redshift 1.5 appears to be the transi- 
tion point between galaxy-magnification dominance and 
magnification-magnification dominance i.e. 2£ g/J > 
for z < 1.5, and vice versa. The inset shows the ra- 



tio £ bs/£gg for several other orientations 9 X (solid lines). 
Comparing the insets for the different redshifts, one can 
see that at low redshifts this ratio quickly drops to unity 
as 6 X moves away from zero, while the drop is much less 
rapid at high redshifts. This is consistent with what we 
have seen in panel (a): the general diffusion of magnifi- 
cation distortion out to larger values of 6x± (and there- 
fore X ) as the redshift increases. The dot-dashed line in 
the inset shows the monopole £ b s divided by £ gg . The 
monopole is defined to be 

,tt/2 

monopole of ^ bs = / £obs(^ x ) sm0 x d# x . (16) 
Jo 
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5x10' 10 2 1.5x10= 

<5x ± [Mpc/h] 

(a) 




90 100 110 120 



<5x [Mpc/h] 

(c) 

FIG. 7: Analog of Fig. [3] for z — 3. The contours in (a) are: 
triple contours), (red) and -0.00002 (green; double contours). 

At a scale of about 100 Mpc/h, the fractional deviation 
of this monopole from £ gs is negligible at z = 0.35, grows 
to a few percent at z = 1.5 and is about 20% at z = 3. 



B. Baryon Acoustic Oscillations 

Panels (c) of Fig. [3] - [7] zoom in on the baryon wig- 
gle for a separation vector oriented along the LOS. Of 
interest, probably in descending order, are the location, 
width and height of the wiggle, and how much they are 
affected by magnification bias. The answers depend on 
whether one looks at £ or, as is commonly done, 6x 2 £. 




o - 

_l i i i i I i i L_ 

50 100 150 

<5x [Mpc/h] 

(b) 




50 100 150 200 

<5x [Mpc/h] 



(d) 

0.01 (magenta), 0.002 (cyan), 0.00045 (black), 0.00015 (blue; 
The dashed line in the lower panel of (d) is £ 99 . 

An important lesson here is that 5x 2 ^ is in fact problem- 
atic: in the presence of magnification bias (and for the 
LOS orientation), the baryon local maximum becomes 
difficult to locate or disappears completely! Confining 
ourselves to £, the fractional changes in the baryon peak 
location, width and height induced by magnification bias 
are shown in the left panel of Fig. [5] for the LOS ori- 
entation, and for several different rcdshifts. Here, we 
go beyond the assumption of (5s — 2)/b = 1 in most of 
the paper, and show these fractional shifts for (5s — 2)/b 
spanning —1 to 2. The shifts can be quite large. They 
can have either sign, depending in part on the sign of 
5s — 2. The magnitude of the shifts depends sensitively 
on (5s — 2) /b and (for width and height) the redshift. All 
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(a)Line-of-sight (LOS) (b)Monopole 

FIG. 8: Fractional shifts in the baryon peak location, width and height as a function of (5s — 2)/b for redshifts spanning 
0.35,1.0,1.5,2.0,3.0. Left panel is for the LOS orientation (where the change in peak location is very similar for different 
redshifts and the redshift labels are therefore omitted), and right panel is for the monopole. The fractional shift is defined 
as (£ bs — £.gg)/£gg for the baryon peak height, and similarly for the peak location and width. The peak width is defined as 
[d 2 (;/dSx 2 /£]~ 1,/2 evaluated at the peak. A note of caution: the precise shifts in these quantities depend on how they are 
extracted from observations. Panels (d) of Fig. [3] - [7] provide a more complete picture of the effects of magnification bias on 
the monopole. 



the shifts vanish at the special point of s = 0.4, where 
magnification bias is absent. One surprise is that even at 
a low redshift of 0.35, where one usually expects lensing 
to be unimportant, the fractional shifts are significant: 
at (5s -2)/b=l, they are - 2%, ~ 50% and - 170% 
for the location, width and height respectively! This is 
fundamentally because the baryon wiggle is a large scale 
feature ~ 100 Mpc/h. At such a large separation, the in- 
trinsic galaxy-galaxy clustering is quite weak. The lens- 
ing correction, on the other hand, is quite large for the 
LOS orientation (9 X = 0) because one galaxy is directly 
behind another. A location shift of ~ 2% is significant 
for baryon oscillation experiments (2l| : it corresponds 
to a ~ 2% shift in the inferred Hubble parameter H(z), 
which corresponds to a shift of ~ 10% for the dark energy 
equation of state w at the redshifts of interest. 

Fig. [5] right panel shows the shifts when one aver- 
ages over directions i.e. the monopole (eq. |16|). The 
shifts are much diminished, but in some cases, still non- 
negligible. For redshifts < 1.5, the location shift is 
< 0.1%, but for redshifts > 1.5, the location shift can 
reach up to ~ 0.6%. The fractional shifts in width and 
height are, as before, a bit larger than that in location. 
An intriguing aspect of Fig. [8]is that, comparing the LOS 
orientation and the monopole, the location shift actually 
differs in sign for (5s — 2)/b > 0. This suggests that the 
monopole is dominated by the transverse modes, which 



is not surprising, given that the monopole downweighs 
the LOS orientation (eq. [To]). 

Note that the monopole can be inferred from data only 
in a cosmology dependent manner i.e. assumptions need 
to be made in relating the observed redshift and angular 
separations to separations in S\ and Sxj_. Ignoring the 
anisotropy induced by magnification bias, it can be shown 
that the baryon peak location in the monopole constrains 
the combination (% 2 /ff) 1 / 3 ([IB]; see Appendix in Pa- 
per II). A ~ 0.6% shift in this quantity translates into 
a ~ 3% shift in the dark energy equation of state w. 
One must keep in mind, however, that in the presence of 
anisotropy, such as that induced by magnification bias, 
(x 2 1 H) 1 ' 3 is not the exact quantity that is constrained 
by the monopole baryon peak. We leave a thorough study 
of this issue for the future. 

One must exercise additional caution in interpreting 
Fig. [5] the precise shifts in the quantities shown de- 
pend on how they are actually extracted from the ob- 
servational data. Typically, one fits some analytic curve 
through the data points to locate and characterize the 
baryon peak. The important point is that magnification 
bias introduces an additive correction to ^ gg which has a 
shape that is uncertain: it depends on b, s and the power 
spectrum. Igoring this correction leads to errors in one's 
fit to the data, and therefore biases in the inferred angu- 
lar diameter distance and Hubble parameter. Precisely 
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how large a bias one would incur is procedure dependent, 
and deserves a careful study. 

Panels (d) of Fig. [3]-[7]emphasize this point by showing 
the difference: monopole of £ b s — £gg , for several different 
values of (5s — 2)/6. (We focus on the monopole instead 
of the LOS orientation, because the monopole is less af- 
fected by magnification bias; the difference £ b s — £, gg is 
quite a bit bigger for the LOS orientation, see panels b 
of Fig. - El) Comparing against the monopole of £ Q b s 
itself, one can see that this difference can be ~ 10% al- 
ready by z ~ 1.5 at the baryon peak. For a limited range 
of scales around the peak, one can perhaps approximate 
this difference by some quadratic function, but the pre- 
cise slope and curvature of this quadratic are uncertain 
and depend on the galaxy/quasar sample in question. 
The issue is whether these should be treated as addi- 
tional parameters in one's fit and how they might affect 
the accuracy of the peak location measurement. We leave 
this for a future paper. 



III. INTERPRETATION 

Fig. [3] - [8] are based on numerically evaluating the 
integrals in eq. (TT2"j) . (TT3|) and (TJjJ) . To gain a deeper un- 
derstanding of what we have seen, it is useful to develop 
order of magnitude estimates for the ratios: i^t/^gg and 



5s 



(1 



</<; 



5s 



50 
1 + 



3 irH A 2 (fc*) 



k* ^^{k^^j 



(17) 



The symbol A 2 (fc) denotes the dimensionless variance at 
scale k and redshift z: Airk 3 P mm (k) / (2ir) 3 . Here, fc** ~ 
1/ y/Sx 2 + 5x\, while fc* is equal to either l/5x±_ or k m , 
whichever is smaller (k m is the scale where k 2 P mm (k) 
peaks; k m > 3 h/Mpc). 

The expressions in eq. (jTTJ) are only approximate. 
Among other things, they fail at separations where 
the correlation functions cross zero or become negative. 
Away from these separations, we have checked that these 
expressions reproduce our numerical results quite well, to 
within factor of a few. Moreover, they make transparent 
several important points. There are two interesting limits 
to consider. 

One is the limit of a small 5\ and a large Sx± (S\ <C 
6x in which case the factors of A 2 's cancel out. It is 
clear that both correlation ratios should then be small, 
as long as 6x± (~ 1/fc*) is a small fraction of the Hub- 
ble scale Hq , and the redshift is not too large. A more 
intuitive way of putting it is that a transverse pair of 
galaxies (i.e. 8\ — and 5x± ^ 0) are on average not 
significantly lensed: for one, they are not lensing each 
other i.e. the galaxy-magnification term vanishes; as far 
as the magnification-magnification term is concerned, it 



is generally small compared to the intrinsic galaxy corre- 
lation unless the redshift is sufficiently high. 

The other interesting limit is that of a large 5\ and a 
small Sx± (5x ^> Sx±), in which case A 2 (fc») ^> A 2 (fc**) 
(because k* ~ 1/Sx± and fc*„ ~ 1/Sx)- This is the limit 
in which the correlation ratios in eq. (TIT)) can potentially 
be of order unity or even larger, with the help of the large 
boost from A 2 (/c*)/A 2 3> 1. What is the physical 
origin of this boost? Consider a pair of galaxies oriented 
along the LOS (i.e. S\ ^ and Sx± = 0). When 5\ is 
large, the intrinsic galaxy correlation is quite weak. The 
magnification bias corrections, on the other hand, can 
be relatively large because this is the case of zero lens- 
ing impact parameter (Sx± — 0). In other words, one is 
fundamentally comparing fluctuations on very different 
scales: large or linear scales in the case of the intrin- 
sic galaxy clustering £ g9 , and small or nonlinear scales in 
the case of the lensing corrections and £ w . The naive 
expectation is that one is performing clustering measure- 
ments on linear scales by considering pairs of galaxies 
separated at a large §x (say 100 Mpc/h), but the truth 
is that the observed clustering is dominated by nonlinear 
fluctuations due to lensing. 

Incidentally, this also means our linear galaxy bias 
assumption is probably not a good approximation for 
pair separations oriented along, or close to, the LOS. 
The galaxy bias is generally expected to first drop be- 
low the linear value for scales just below the nonlin- 
ear scale, and then climb when approaching zero lag 
[TTI | . This means the effects of gravitational lensing on 
the LOS orientation should be enhanced by a nonlinear 
galaxy bias compared to our predictions. One can see 
this by recalling the expression for the observed correla- 
tion: £ bs = £ gg + 2£ 5Al + In the LOS orientation, 
and for a large separation, the relevant bias for £ g9 is the 
linear galaxy bias b that we have been using all along, 
while the relevant bias for (at Sx± = 0) should re- 
ally be a nonlinear and presumably enhanced galaxy bias 
> b. Therefore, the effects of gravitational lensing are in 
fact underestimated by our calculations which use a lin- 
ear galaxy bias. The precise scale dependence though 
of the nonlinear galaxy bias is quite complicated and is 
highly sample dependent. We hope to pursue a thorough 
study of the effects of nonlinear galaxy bias in the future. 
It is important to stress that the distinctive anisotropy 
displayed in eq. (|15p holds true even in the presence of 
nonlinear galaxy bias - this is because eq. (fT5")) follows 
from eq. (HU), (flU)) and (fLT| which make no assumptions 
about galaxy biasing. 

Irrespective of details of the nonlinear galaxy bias, the 
following statement is expected to hold: magnification- 
bias has the strongest effect on the observed correlation 
function when the separation (Sx has a large magnitude 
and points close to the line- of- sight direction. 

Eq. (fT7|) also tells us that £ MAt is generally expected to 
be larger than £ gAt when the redshift is sufficiently large. 
Sufficiently large here means when xHq > 1. This occurs 
at z ~ 1.5 for the cosmological parameters we adopt. 
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Conversely, is generally quite small at low redshifts 
because of the cubic power in the factor (x^o) 3 - I n this 
case, £ gfJ , is more favored if the LOS separation 5\ is large 
enough. 



IV. DISCUSSION 

Our findings can be summarized as an opportunity and 
a challenge. 

Opportunity. Gravitational lensing, through magnifi- 
cation bias, introduces a distinctive anisotropy to the ob- 
served 3D galaxy/quasar correlation function (eq. |15|). 
We have studied its shape in £111 Al The correlation 
is preferentially enhanced (in the positive direction for 
s > 0.4) in the LOS orientation, vaguely resembling the 
finger-of-god (FOG) effect due to virializcd peculiar mo- 
tions. However, as will be discussed in detail in Paper II, 
the precise shape of magnification distortion differs from 
that of FOG or redshift distortion in general - redshift 
distortion does not have the distinctive linear dependence 
on the LOS separation that £ gM has. The distinctive 
lensing induced pattern creates an interesting opportu- 
nity: it is in principle possible to separately obtain from 
data all three contributions: the galaxy-galaxy, galaxy- 
magnification and magnification-magnification correla- 
tions (Fig. [2]). This generalizes earlier work on mag- 
nification bias [1, 0, !, H El which focused on angular 
correlations in the large redshift separation limit, where 
only the galaxy-magnification correlation survives. Fu- 
ture galaxy/quasar (spectroscopic or photometric) red- 
shift surveys can be used to measure the galaxy-galaxy, 
galaxy-mass and mass-mass power spectra, even without 
measurements of the galaxy shapes. 

Challenge. The same effect must be accounted for in 
interpreting galaxy/quasar clustering data. This presents 
an interesting challenge to precision measurements, such 
as those that hope to use the baryon oscillation scale to 
yield stringent constraints on dark energy. Contrary to 
naive expectations, magnification bias is not necessarily 
negligible at low redshifts. In the LOS orientation, where 
its effects are largest, the shift in the baryon oscillation 
scale is remarkably insensitive to redshifts, from z ~ 0.35 
to z ~ 3 (Fig. [8] panel a). The shift could reach up to 
~ 2 — 3%, which translates to a ~ 10 — 15% shift in 
the dark energy equation of state from the inferred Hub- 
ble par ameter. Existing baryon oscillation measurements 
[l5l. Till [ttI. [l8l IT9II have not reached this kind of accuracy, 
but the future ones will [2l|. The shift for other orien- 
tations are smaller. We have considered the shift in the 
monopole and find that it is much diminished, reaching 
up to ~ 0.6% at z ~ 3 (Fig. [8] panel b). However, one 
should keep in mind that the precise shifts really depend 
on exactly how the baryon oscillation scale is extracted 
from data. The point is that magnification bias intro- 
duces corrections to the correlation function that are de- 
pendent on scale, bias and the number count slope, for all 
orientations and for the monopole (Fig. [3] - [7] panels d). 



The question is when data are fitted without accounting 
for these corrections, what kind of bias would one incur 
in the inferred parameters? 

This leads us naturally to several interesting questions 
to be explored in the future. 

On the opportunity side, what is the optimal scheme 
to extract the galaxy-galaxy, galaxy-magnification and 
magnification-magnification correlations from realistic 
data? Fig. [2] represents a proof of concept but is likely 
not the optimal method. How to best weigh the relative 
contributions of different galaxies with different bias and 
number count slope? Also, how to best guard against 
the possibility of confusion with dust extinction? The 
distinctive dependence of the magnification corrections 
on the number count slope can probably be exploited for 
this purpose [8j. 

On the challenge side, how would magnification bias 
impact baryon oscillation measurements if one accounts 
for precisely how the oscillation scale might be extracted 
from future data? How might the impact be different 
in Fourier space versus real space? We will address this 
question partially in Paper II. 

Our calculations assumed a constant bias b which is 
certainly incorrect at some level. A complete calculation 
that accounts for nonlinear bias and distortions of all 
types (peculiar motions, Alcock-Paczynski and magnifi- 
cation bias) should be made to properly interpet preci- 
sion galaxy clustering measurements [29| . Lastly, in light 
of our findings, it is probably interesting to revisit some 
longstanding puzzles, such as the well known results of 
the pencil beam surveys by [3(| HH . The enhancement of 
the LOS correlation by magnification bias offers a plau- 
sible explanation for the excess correlations, though not 
the periodicity, seen in [3(J. With the magnitude cut-off 
used by [3(j, it is likely that (5s — 2)/b ~ 1 (see also 
[32l]). in which case we can scale the result of Fig. [3lc) 
to the mean redshift of ~ 0.2 using eq. (jTTJ) and obtain 
an enhancement factor of ~ 2 at the scale ~ 100 Mpc/h. 
Coupled this with the fact that the baryon wiggle was 
not taken into consideration back in the days of [30| , the 
total enhancement factor by magnification bias + baryon 
oscillations is ~ 4. Moreover, as noted in £11111 accounting 
for nonlinear galaxy bias in the lensing correction would 
likely lead to a further enhancement factor of a few. The 
net enhancement factor could well be ~ 10, making the 
excess correlations seen in [30] within reach. It would 
be worthwile to carry out a more careful analysis, taking 
into account properties of the precise galaxy sample in 
question, and the actual widths of the pencil beams. It 
might also be interesting to revisit some of the measure- 
ments of Lyman-break galaxy clustering at z ~ 3 [33| in 
the light of our findings. 
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APPENDIX A: MAGNIFICATION BIAS AND 
OTHER LENSING CORRECTIONS TO THE 
OBSERVED GALAXY CLUSTERING 

In this Appendix, we attempt to be as general as possi- 
ble in considering the effect of gravitational lensing on the 
apparent spatial clustering of galaxies. In the process, we 
will rederive the magnification bias effect, generalizing it 
to the case of a gradual rather than sharp magnitude cut- 
off. We will also discuss the role of stochastic deflections 
which contribute to higher order corrections. 

Let z, 0g)df g be the intrinsic number density of 

galaxies at redshift z, source (unlensed) position 9 g , and 
an unlensed flux in the range f g ± df g /2. We use the 
subscript g to denote quantities that are intrinsic to the 
galaxies, or in other words, prelensed/unlensed quanti- 
ties. Note also that the number density here can refer 
to either an angular number density or a volume number 
density. Gravitational lensing introduces the transfor- 



mations f g — > /, 9g 
following relations: 



9 and $ 9 — > $, which satisfy the 



where 



<S>(f,z,9)dfd 2 9 = <S> g (f g ,z,9 g )df g d 2 



e„ = e + se 



(Al) 



(A2) 



dct. 



f = Af 9 
= 1/ A 



The quantity A denotes the magnification in the observed 
direction 9, and 69 is the lensing displacement. 

A given galaxy sample can be modeled by an efficiency 
function e(/) such that the observed galaxy density is 



n(z,9) = J e(f)$(f,z,0)df 



(A3) 



The simplest example of e(f) is a step function which 
equals unity if / > / m in. and vanishes otherwise: 



<f) = e(/ min . < /) 

but let us keep e(f) general. 

Using eq. (|Aip and (IA2|) . we can see that 



(A4) 



n(z, 0) = J e(Af g )$ g (f g , z,9 + S9)df g (A5) 



Note that the lensing magnification A is a function of 
both the source redshift z and the direction 9. 

Our expressions so far are completely general: no weak 
lensing approximation has been made (aside from not 
summing over multiple 69 's to account for the possibility 
of multiple images). One could if one wishes use eq. (|A5j) 
as the starting point for investigating the lensed galaxy 
clustering. The discussion is much simplified, however, 
under the weak lensing approximation i.e. A ~ 1 + 2k, 
where \k\ -C 1. In that case, 



where 



n(z, 9) = n g (z, 9 + 89) [1 + (5s - 2)k] (A6) 



<f 9 )%(f 9 )df g (A7) 



is the galaxy density if lensing magnification were absent 
(we have suppressed the z and 9 + 59 dependence), and 



[2.5U,]- 1 J ^f g ^ g (f g )df g 



(A8) 



The definition of s might look a bit unfamiliar, but if 
e does indeed take the form of the step function (eq. 
|A4| ) , s is the (magnitude) slope of the cumulative num- 
ber counts at the faint end cut-off: 



dlogi 



dm 



(A9) 



where m = — 2.51og 10 / and the derivative is evaluated 
at m= -2.51og 10 / m i n .. 

As long as (A) = 1 (or (k) — 0) which holds if multiple 
imaging can be ignored, (n) = (fi g ). Let us call this mean 
density n. The fluctuation in the observed galaxy density 
8 = n/n — 1 is related to the fluctuation in the intrinsic 
galaxy density 8 g = n g /n — 1 by 



8(z,0) = 8 g (z,d + 80) + 

(5s - 2)k(z, 6) x (1 + 8 g (z, 9 + 89)) 



(A10) 



This is equivalent to eq. ^ and used in the rest 
of the paper, except that we have consistently ignored 
the effects of stochastic deflection 89, and we have con- 
sistently neglected nS g as small. In other words, gravi- 
tational lensing strictly speaking should have 2 qualita- 
tively different effects on the observed galaxy clustering. 
Magnification bias is embodied in the term (5s — 2)k, 
which accounts for the effect of an overall magnification 
(or demagnification) on the number counts. Stochastic 
deflections, on the other hand, introduces a remapping 
of the galaxy density field. The remapped field can be 
Taylor expanded as 



8 g (9 + 89) ~ 8 g {9) 



(All) 



The approximation we have been making is in effect a 
linear approximation: the second term on the right is 
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second order in perturbations and is therefore ignored 
(just as we throw away the n8 g term also). Similarly, the 
s defined in eq. (|A8j) can be replaced by one with f g —*f 
and $ 9 —> <I>. Corrections to this when multiplied by re 
as in eq. (|A10[) are second order. 

It is worth emphasizing that the gravitational lens- 
ing effect on diffuse backgrounds, such as the microwave 
background or redshifted 21cm background, is quite dif- 
ferent [H, [H[ . There, because the observable is specific 
intensity which is conserved by lensing, the only lensing 
effect is remapping by stochastic deflections i.e. there is 
no analog of the magnification bias term (5s — 2)re. 

One might worry that in the nonlinear regime where 8 g 
is not small, one should not ignore the effects of stochastic 
deflection relative to magnification bias. This certainly 
deserves more study. Zahn & Zaldarriaga [HJ carried out 
a detailed study of the impact of stochastic deflection on 
the redshifted 21cm background, including all higher or- 
der terms. Stochastic deflection appears to have a rather 
small overall impact, suggesting that ignoring it might 
not be a bad approximation. 



APPENDIX B: HIGHER ORDER CORRECTIONS 

Eq. (TT2")) . (fP3")) and [pH]). and the corresponding expres- 
sions in Fourier space, are written down by keeping only 
the dominant terms in an expansion using (xi — x)lx or 
(X2 — X)/x as a small parameter. It is worth discussing 
how large the corrections are expected to be. One might 
naively expect the corrections to be one order higher: for 
instance, that the dominant correction to £^(1 ! 2) should 
be of the order of what is given in eq. (|13p multiplied by 
(Xi — X2)/x- This expectation turns out to be false for 



£ MM and £ gg , but valid for Let us see how this comes 
about in the case of Taylor expanding eq. flT]) for 
£ MM (1;2) in small 5xi = Xi - X and 6x2 = X2 - X, w e 
have 



U(Xi,0i;Xa,0a) = [^n m (5s - 2)] 1 



(Bl) 



d X ' 



(x - x')x' 



f (Pkj_ p , , , s ikL-x' (01-02) 

J ^27r) 2 lm 

[l+0(5xi/x) + 0(S X 2/x)} 



where the terms 0{8 X i/x) an d 0(6x2/ X) signify all the 
next-to-dominant order corrections which in general do 
not vanish. The important realization is that the term 
0(8xi lx) an d t ne term 0(8x2 /x) have exactly the same 
coefficients, and so they sum to zero as long as 8x1 + 
8x2 = i.e. x is the mean of X i an d X2- Therefore, one 
has to go one order higher in the Taylor expansion to fig- 
ure out the actual corrections. A similar argument works 
for the galaxy auto-correlation t; gg . The magnification- 
galaxy cross-correlation £ SAt (1; 2) +£ gfl (2; 1), on the other 
hand, turns out to be different: in this case, the anal- 
ogous terms 0(8xi/x) an( i 0(8x2/ X) have different co- 
efficients. In general, when the pair (1;2) is averaged 
over some survey volume, one expects the (fractional) 
correction to our expression for the magnification-galaxy 
cross-correlation to be 0(\xi — X2|/x)- If one is worried 
about this correction, there is no reason why one cannot 
use instead the exact expression (within Limber approx- 
imation) for £ gAt in eq. ©. 
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